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ABSTRACT 

In wind-fed X-ray binaries the accreting matter is Compton cooled and falls freely onto the 
compact object. The matter has a modest angular momentum I and accretion is quasi-spherical 
at large distances from the compact object. Initially small non-radial velocities grow in the 
converging supersonic flow and become substantial in the vicinity of the accretor. The stream- 
lines with I > (GM i?*) 1 / 2 (where M and are the mass and radius of the compact object) 
intersect outside and form a two-dimensional caustic which emits X-rays. The stream- 
lines with low angular momentum, I < (GMR*) 1 / 2 , run into the accretor. If the accretor 
is a neutron star, a large X-ray luminosity results. We show that the distribution of accretion 
rate/luminosity over the star surface is sensitive to the angular momentum distribution of the 
accreting matter. The apparent luminosity depends on the side from which the star is observed 
and can change periodically with the orbital phase of the binary. The accretor then appears as 
a 'Moon-like' X-ray source. 
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1 INTRODUCTION 

Wind-fed accretion is believed to occur in massive X-ray binaries 
(see e.g. King 1995 for a review). The massive donor companion 
(OB star) produces a substantial wind, up to M w ~ lO _6 M0/yr, 
which is partly captured by the compact companion. The wind ma- 
terial is captured from an accretion cylinder of radius 



Ra — 



2GM 



(1) 



where w is the wind velocity and M is the mass of the accre- 
tor (Hoyle & Lyttleton 1939, Bondi & Hoyle 1944). The typ- 
ical w ~ 10 8 cm s _1 for OB stars, so that R a ~ 3 x 
1O 1o (M/M )(w/1O 8 )" 2 cm. If the wind is isotropic then the ac- 
cretion rate is M ~ (l/4)(R a / A) 2 M w where A is the binary sep- 
aration. The accretion rate is substantial in close binaries only, with 
orbital periods P ~ a few days. The captured fraction, M/M m , 
can be increased if the donor is a Be star which has a prominent 
slow equatorial wind (see e.g. van Paradijs & McClintock 1995) or 
if the wind is prefocused by the tidal effects (Blondin, Stevens & 
Kallman 1991). 

Gas captured from the accretion cylinder falls many decades 
in radius down to the radius of the compact object R, where 
the X-rays are produced. R* equals r g — 2GM/c 2 ~ 3 x 
10 5 (M /Mo) cm for a black hole (BH) and about 3r 9 for a neu- 
tron star (NS). Owing to orbital rotation of the binary, the captured 
gas possesses a net angular momentum with respect to the accretor. 
The average angular momentum 1 can be estimated (Illarionov & 
Sunyaev 1975; Shapiro & Lightman 1976). It is directed perpen- 



dicularly to the binary plane and equals l z = ("(1 / '4)fi7?2 where 
Q. — 2n/P is the angular velocity of the binary and the numerical 
factor ( 1 depends on adopted assumptions (see e.g. Wang 1981; 
Livio et al. 1986; Ruffert 1997, 1999). The angular momentum is 
small and the infall is radial at R 3> R* ■ 

Compton cooling by the central X-ray source makes the inflow 
highly super-sonic inside the Compton radius, Rc ~ 10 10 cm < 
R a (see Illarionov & Kompaneets 1990 and Section 2.2). Initially 
small non-radial velocities v±_ — l/R grow in the freely falling 
flow and exceed the radial velocity component v r ~ (GM / R) 1 ' 2 
at Rd ~ I 2 /GM. Accretion can be assumed to be radial if Rd <C 
R,, i.e. if I <C I* where 



I, = (GMR, 



,1/2 



(2) 



The deviations from the radial pattern are important if I is compa- 
rable to h. 

In BH and NS binaries, Z» ~ r g c and l z /r g c w 
1.5C-P _1 (A//M Q )('u;/10 8 )" 4 where the binary period P is mea- 
sured in days. The observed P in massive (OB) X-ray binaries is 
typically a few days (White, Nagase & Parmar 1995; Tanaka & 
Lewin 1995). One thus concludes that l z is about Z» in these sys- 
tems. Note that the angular momentum of a particular streamline 
varies substantially around the average value. For instance, if the 
flow is in solid body rotation at R i?* then I is highest for 
streamlines in the equatorial plane and vanishes on the polar axis. 

Matter with I < U runs directly into the accretor before it 
reaches the equatorial plane. If the accretor is a neutron star then 
a strong shock results and X-rays are produced (Zel'dovich & 
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Shakura 1969; Shapiro & Salpeter 1975). The surface brightness of 
the star is determined by the distribution of the accretion rate over 
its surface, AM /AS. In this paper, we find that AM /AS is sensi- 
tive to the angular momentum distribution in the flow. We assume 
a weakly magnetised NS (B & 10 8 G), so that the magnetic field 
does not affect the ballistic trajectories of the freely falling matter. 
The resulting surface brightness of the star is inhomogeneous and 
the apparent luminosity depends on the side from which the star is 
observed. The apparent luminosity can then change as the binary 
executes its orbital period (we dub it the 'Moon' effect). 

By contrast, if the accretor is a black hole then matter with 
I < U plunges into the event horizon without producing substantial 
emission. 

The streamlines with I > Z* intersect in the equatorial plane 
(the plane of symmetry) at R > R, . The loci of the intersections 
form a two-dimensional caustic. If the accretor is a black hole then 
the caustic is the only source of X-rays from the accretion flow. 

The paper is organised as follows. In Section 2 we briefly re- 
view the pattern of wind-fed accretion on large scales, at distances 
~ R a from the accretor. In Section 3 we write down the equations 
of the freely falling flow inside the Compton radius. In Sections 4 
and 5 we focus on the very vicinity of the compact object. We dis- 
cuss asymmetric accretion onto the surface of a NS (Section 4) and 
then caustics outside the accretor (Section 5). 



2 WIND-FED ACCRETION ON LARGE SCALES 
2.1 The trapping of the wind matter 

The radius of the accretion cylinder is small compared to the binary 
separation, R a 5s 10 11 cm <C A ~ 10 12 cm, and hence the flow 
in the cylinder is nearly plane-parallel before it gets trapped by the 
gravitational field of the accretor. The wind reaches the compact 
companion on a time-scale A/w ~ 10 4 s which is much shorter 
than the orbital period P ~ 10 6 s. As a first approximation, one can 
assume the accretor to be at rest and the flow to be axisymmetric 
around the line connecting the two companions. 

Let us introduce coordinates (a;, y, z) so that the y-axis is di- 
rected from the donor to the accretor and the z-axis is perpendicular 
to the binary plane, and choose the coordinate origin at the location 
of the accretor (see Fig. 1). Each streamline of the flow is specified 
by two impact parameters (xo, zo) at R 3> R a . The initial angular 
momentum of a streamline (xo, zo) is 1 = (— zqw, 0, xow) and its 
absolute value equals Z = bw where b — (a; 2 , + Zq) . The net 1 
integrated over a ring (b, b + db) vanishes, which is a consequence 
of the assumed symmetry around the y-axis. 

For an initially super-sonic wind, a bow shock forms at dis- 
tance ~ R a from the accretor. Hunt (1971) first studied gas dy- 
namics behind the shock and showed that a spherically symmetric 
inflow forms at R < R a (see also Petrich et al. 1989; Ruffert 1997, 
1999). Fig. 1 shows the picture of accretion. The transformation 
of the uniform plane-parallel flow into the isotropic spherical infall 
can be described as follows. A streamline with an initial impact pa- 
rameter b eventually infalls radially at some angle 3 with respect to 
the — y axis, 

/ b \ 2 1 — cos 3 ,o 2 2 

This transformation induces a map (xo, zo) — * (0, ¥>), 
Xo sin 6 sin ip z cos 6 

Ra ^/2(1 + sin 9 cos ip) ' R a ^2(1 + sin cos 




y 



Figure 1. Schematic picture of wind-fed accretion. The initially plane- 
parallel flow is captured by the accretor and transformed into a spherical 
inflow. Soo denotes a sphere of radius R ~ Rq (eq. 7). Inside Sao, the 
flow is falling freely (Section 3). The coordinate origin is at the location of 
the accretor, A. 

Here 6 is the polar angle measured from the z-axis and ip is the az- 
imuthal angle measured in the xy plane from the (— jy)-axis. Note 
that the boundary of the accretion cylinder x\ + Zq — R a trans- 
forms into one point 6 = 7r/2, ip — n. 

The mapping (3) assumes that the flow is laminar behind the 
bow shock. Numerical simulations show that the flow is unstable 
if the bow shock is strong, with a high Mach number. However, in 
the case of modest Mach numbers, the fluctuations are weak and 
the flow is approximately laminar (e.g. Blondin et al. 1990; Ruffert 
1997, 1999). This is the most likely case if accretion occurs in the 
radiation field of a luminous X-ray source (see below). Note also 
that the streamlines cross the shock nearly normally (Hunt 1971); 
therefore the shock does not generates vorticity in the flow (Landau 
& Lifshitz 1987). 

2.2 Compton heating/cooling 

If the central X-ray source is luminous, L > 10~ 3 Le where 
Le = 4ncGMm p /<jT is the Eddington luminosity, then Compton 
heating/cooling affects strongly the dynamics of accretion (Ostriker 
et al. 1976; Illarionov & Kompaneets 1990). We here discuss two 
effects: (i) On large scales, R J; R a , the wind matter is preheated 
and its Mach number is reduced to M ~ 1. (ii) On small scales, 
R < Rc (see eq. 7), the accreting matter is cooled by the X-rays 
and falls super-sonically (Zel'dovich & Shakura 1969). 

2.2.1 Preheating of the wind 

The initial temperature of the wind is To ~ 10 5 K and the typical 
Mach number is Ma = w/c s ~ 10 where c 3 = (WkT/Sirip) 1 ^ 2 
is the sound speed. One can evaluate how M decreases along 
a streamline when approaching R ~ R a . The Compton heat- 
ing/cooling is dominant in the energy balance of the highly ionised 
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(5) 



plasma in the accretion cylinder (see e.g. Igumenshchev, Illarionov 
& Kompaneets 1993). Then the energy balance reads 

where T and v are the plasma temperature and velocity, respec- 
tively, Tc ~ 10 8 K is the Compton temperature of the X-ray 
source, tc = 3irm e c 2 R 2 far L is the time-scale for Compton cool- 
ing, and L is the source luminosity. Since Tq 3> T at R R a , the 
Compton heating term is dominant and we keep only this term on 
the right-hand side of the energy equation, 

AT _ Tc 
At ~ t c ' 

In the upstream region, the wind matter is falling freely. From 
the angular momentum conservation we have 

R 2 ^=bw 
at 

where b is the impact parameter of the streamline and a is the an- 
gle between the radius vector R and the (—y) axis (see Fig. 1). 
Combining this equation with the energy equation (5), we get 

AT_TcR 2 __ a T LT c 
Aa tc bw 3irm e c 2 bw 

Taking R — bj sin a w b/a, we get the temperature 

T*T(R)= aTLT C . 

3nm e c 2 Rw 

At R = R a it gives a Mach number 

2 _ 3m p w 2 _ w f L \~ l (T c \~ l 
VO.ILeJ U0 8 J 



const. 



(6) 



M 



lOfcT 



10 s 



As a result of the Compton preheating, M decreases markedly be- 
low the initial Mo, and, correspondingly, the strength of the bow 
shock is reduced. In bright hard sources, the heated wind may be- 
come subsonic atR^R a and then the shock disappears. 

2.2.2 Compton radius 

Compton heating leads to an inflow-outflow pattern of accretion 
with M. ~ 1 down to the Compton radius (Illarionov & Kompa- 
neets 1990; Igumenshchev, Illarionov, & Kompaneets 1993), 



GMm p 



3 x 10 J 



M 



(7) 



Typically, Rc is ~ 10 times smaller than R a . 

Inside Rc, a spherical inflow forms. Its temperature exceeds 
Tc , so that the gas is cooled by the X-rays rather than heated (see 
eq. 4). The cooling leads to a high Mach number of the spherical 
inflow, M~>1. 



2.3 The trapped angular momentum 

When the accretor orbital motion is taken into account, the flow 
pattern is no longer symmetric around the y— axis and the captured 
matter should have a small net angular momentum l z ~ QR 2 with 
respect to the accretor. One would like to know the distribution of 
I around l z in the accretion flow. This distribution governs the flow 
dynamics in the vicinity of the accretor where deviations from the 
spherical pattern become substantial. 

Assuming a weak bow shock, the flow is almost laminar all 
the way. Then the transformation of the accretion cylinder into 



the spherical inflow is given by equation (3). The corresponding 
transformation of angular momentum is not known. The analysis 
of small perturbations in a converging flow (e.g. Lai & Goldre- 
ich 2000 and references therein) shows that the rotational mode, 
«rot oc r _1 , has the fastest growth inwards, while sonic modes are 
damped. The rotational mode is probably dominant on Rc and we 
restrict our consideration to this mode and the associated angular 
momentum. 

A particular streamline with impact parameters (xo, zo) con- 
tributes to the net trapped angular momentum proportionally to 
Six 2 , (see e.g. Shapiro & Lightman 1976). If this magnitude con- 
serves along a streamline then mapping (3) also determines the dis- 
tribution of the trapped I over angles 9, ip at R <C -Rc, 



l z {6, <p)=lo 



sin 2 9 sin 2 ip 



I - l ± 



(8) 



1 + sin 9 cos (fi 

We assume that the infall angular momentum is associated with 
rotation around the 2-axis only, i.e. we assume v v » vg. Then the 
orbital angular momentum of a streamline is 



I = 



(9) 



The corresponding non-radial velocity is v± = v,p — l z /R. The 
distribution (8) is not axisymmetric, which leads to an essentially 
three-dimensional pattern of accretion in the vicinity of the com- 
pact object. 

In the other limiting case, the trapped / is efficiently redis- 
tributed between the streamlines so that they come to solid body 
rotation around the z-axis with a common angular velocity u. Then 



l z (0) = l sm 2 e, 



I = flo. 



(10) 



The non-radial velocity is v± — v v — luR sin 9 = l z /R and the 
infall is axisymmetric. 



3 FREE FALL 

In this section, we write down the equations describing the matter 
free fall in Newtonian gravity (the relativistic equations are given in 
an accompanying paper, Beloborodov & Illarionov 2000). We con- 
sider sub-Eddington sources only, where radiation pressure does 
not affect the flow dynamics. We are interested in the part of the 
ballistic trajectories before the collision with the accretor or the 
intersection with the equatorial plane, the symmetry plane of the 
flow. 

At R <§C Rc, the free fall is nearly parabolic. A parabolic 
trajectory is described by the equation 



l 2 R* 



h = Vgmr, 



(11) 



H(l — costp) ' 

where 1 is the orbital angular momentum of a streamline and ip is 
the angle between the changing radius-vector of the streamline R 
and its initial radius- vector at infinity Roo (see Fig. 2). The label 
'oo' corresponds to distances R ~ Rc. 

The angular distribution of the accretion rate at a sphere Sr of 
radius R is different from the initial uniform distribution at Soo ■ the 
rotation defocuses the inflow and makes it non-uniform. A stream- 
line that starts at Oao, ipoo will cross Sr at 0, ip which satisfy the 
relations (see Fig. 2) 



COS 9 = COS 9oo cos ip, 

. , ., simp 

sm(<p - Lpoc) = — -, 



(12) 
(13) 
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4 NEUTRON STAR AS A MOON-LIKE X-RAY SOURCE 




Figure 2. Geometry of the free fall. At point A (radius r = (l/l,) 2 R t ), 
the ballistic trajectory intersects with the symmetric trajectory from below 
the equatorial plane. If r < then the trajectory instead runs into the 
accretor, before it can reach the equatorial plane. 

where cosi/> — 1 — l 2 R*/l 2 R (from eq. 11). The accretion rate 
distribution at Sr is determined by the Jacobian of the mapping 
(#oo, <Poo) — » (0, <p). After some algebra we get the Jacobian 



A _ d{co S 6,<p) 

d(cos 8oC ; (PoO 

where 

dcostp 



cos V>— ctg# c 



d cos V> ctgtp d cos ip 
ouoo sin Uoo dip 



(14) 



d0 o 



21R* dl 
l*R dOoc 



dcostp 



2lR t dl 
' HR 



□ ' 9v?oo UR d<poc ' 

The mapping is one-to-one if A > at any 9^ , Lp x . Vanishing 
of the Jacobian implies intersection of the ballistic trajectories. The 
streamlines then approach each other and the pressure effects must 
switch on and prevent the streamlines from intersecting. If the ve- 
locity of the approaching is super-sonic then shocks must occur. 
It is convenient to rewrite equation (14) as 



A 2 ) ll + q 



d\ 

cVoc 



+ Ctgflc 



d\ 2 

'<96>oo 



(15) 



where A = {l/l„){R„/Rf/ 2 and q = 2(2 - A 2 )~ 1/2 sin -1 8^. 
In the axisymmetric case (Z does not depend on ipoo) the Jaco- 
bian is positive if Al / Asm8 x > 0, i.e. if 1(8 oc) increases to- 
wards the equatorial plane. Then the ballistic approximation can 
be used down to the surface of the accretor (in the case I < Z*) or 
down to the equatorial plane (in the case Z > Z„). The condition 
dZ/d sin 6>oo > is satisfied for e.g. Z-distribution (10). 

With A > we have a simple expression for the angular dis- 
tribution of the accretion rate M over a sphere of radius R 



AM 
AO. 



1 AM 
A dfioo 



AO, = dcos 9A(p. 



(16) 



In this section, we study inflows with sufficiently small angular mo- 
mentum, Z < Z*, which have no caustics outside the accretor. Each 
ballistic streamline runs into the accretor. Throughout this section 
we will assume that the accretor is a (weakly magnetised) neutron 
star. The collision of accreting matter with the star is accompanied 
by a strong shock. We assume that the shock is radiatively efficient 
and that it is held down to the star, i.e. the height of the shock is 
small compared to i?* . This situation is likely to take place at suf- 
ficiently high accretion rates (see Shapiro & Salpeter 1975). Then 
the accreting matter is in free fall until it reaches the surface of the 
star. 



4.1 Surface distribution of M 

The presence of non-zero angular momentum leads to an inho- 
mogeneous distribution of M over the surface of the star. We 
now compute this distribution. With a homogeneous accretion rate 
through Soc, AM /AQ.OO — M t ot/47r, equation (16) yields 



AM ia \ 



M u 



4tyR 2 A' 



(17) 



First consider the axisymmetric case and for illustration take 
the Z-distribution (10) which corresponds to solid body rotation at 
R — ► co. The accretion rate AM /AS is then given by 

AS U ~ 47ri? 2 [l + (Z /Z») 2 (3cos 2 (9 



1)]' 



The flux of matter impinging the accretor increases towards the 
equatorial plane (Fig. 3a). This is a consequence of the fact that 
sin 8 increases along a rotating trajectory, so that the flow gets con- 
centrated towards 8 = n/2. In the example shown in Fig. 3a the 
accretion rate at the equator is twice as large as that at the polar 
cap. 

In general the inflow does not need to be axisymmetric. In 
particular, the trapped Z in the 'tail' (<Poo ~ 71") may differ from 
Z trapped at ipoc ~ 0. For example, the distribution (8) gives the 
whole set of <£>-harmonics, m = 1, 2, .... To study the effects of the 
(^-asymmetry, we consider here the first (m = 1) mode, i.e. 



l(8 x ,<p a 



) sinS 00 (l — acos<y3oc) 



a < 1. 



(18) 



The condition that Z < U over Soo requires Z (l + a) < Z„. Under 
this condition, equation (15) yields A > at R > _R*, i.e. there 
are no intersections of the ballistic trajectories outside the accretor. 

We have computed AM /AS on the surface of the star for in- 
flows with Zo = 0.5Z* and three values of a = 0, 0.5, and 1 (Fig. 3). 
At a ^ the accretion rate is no longer axisymmetric. The inho- 
mogeneity of AM /AS increases with increasing a. One can see two 
effects: 

(i) The accretion flow concentrates towards the equatorial plane 
along the streamlines with the highest angular momentum (the 
maximum Z = Z max is at 8^ = 7r/2 and ifoo = ir). For Zo = 0.5Z* 
and a = 1 we have Z max = U and the corresponding streamline 
marginally touches the star at 8 = 7r/2 and ip — 3n/2. Here A = 
and dM/dS — » oo (Fig. 3cd). This critical point is a 'seed' of the 
outward caustic: further increase in Zo would result in Z max > Z* 
and then the streamlines collide in the equatorial plane outside the 
star (see Section 5 and Fig. 5). 

(ii) The flow concentrates around the meridian ip ~ 77r/4 pro- 
ducing a 'Moon-like' spot on the accretor. The spot is formed by 
the trajectories starting at tpac ~ 37r/2 where dl/dtpoo reaches 
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Figure 3. The distribution of dM/d,S over the surface of the accretor. 
Dark regions correspond to high dM /d,S. Plots (a-c) show the view from 
the top (the accretion flow is symmetric about the equatorial plane). The 
angular momentum of the flow is given by equation (18) with lo = 
0.5(GMi?*) 1 /2. pi ots ( a ) > (b), and (c) show the cases a = 0, a = 0.5, and 
a = 1, respectively; (d) is same as (c) but the moving star is viewed from 
the head-on direction in the equatorial plane. Solid curves are the contours 
of constant dM/dS with a logarithmic step 0.15 in plots (a) and (b), and 
0.3 in plots (c) and (d). The label marks contours dM /d,S = Mtot/47r 
(A = 1). The arrows show the projections of free-fall streamlines onto 
the star surface for 8^ = 7r/6, it/3 and ipoo = 0, tt/2, it, 3tt/2 (see 
eq. 12,13). 

its minimum — a/osinf?oo (see eq. 18). Angular momentum de- 
termines the velocity of rotation in the p— direction. Like usual 
one-dimensional motion with position-dependent velocity, the flow 
with dl/difoo > diverges in the (^—direction {dM/dip de- 
creases) while the flow with dl/difoo < converges {dM/dip 
decreases). Therefore dM/dS peaks at the trajectories with mini- 
mum dl / dip oo ■ These trajectories start at ipoa — 3tt/2 with angular 
momentum lo sin 0^, and they collide with the star at ~ 77r/4. 

The Moon effect is clearly seen from a simple analytical con- 
sideration. Assuming a small lo and neglecting the A 2 terms in 
equations (13,15) one gets 

p = pec + y~V2(1 — acos<,0oo) , 

A = 1 + a^-v^sin^oo. 

The dependences ip(ipoo) and A(<^oo) describe a cycloid A(tp) 
in its standard parametric form. The peak of the accretion rate 
dM/dS = (M t ot/47T^)[l - a(7 /'*)V2] _1 is achieved at the 
meridian ip = 37r/2 + (£n/Z*)v2. In the case lo = 0.5Z* it yields 
p = 270° +40°. 

4.2 The observed luminosity 

On the NS surface, the flux of accreting mass is converted into ra- 
diation with efficiency GM/ R,c 2 . The emerging radiation flux is 
given by 
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Figure 4. The apparent luminosity L b s normalised by the total true lumi- 
nosity L = GMMtot /R* ■ (a) - i bs as a function of the binary inclina- 
tion i in the case of a = 0, lo = 0.5£„ (see eq. 18). In this case, L b s does 
not depend on the orbital phase of the binary <j> owing to the axial symmetry 
of the accretion flow, (b) - L ohs as a function of <j> in the case of a = 0.5, 
Iq = 0.5i*. The solid, long-dashed, short-dashed, and dotted curves corre- 
spond to inclinations i = 7r/2, tt/3, it/6, and 0, respectively, (c) - Same 
as (b) but with a = 1. 



F(e,<p) 



GM dM 
~R7~dS' 



(19) 



The distribution dM /dS shown in Fig. 3 also displays the bright- 
ness distribution as seen by an observer located at the polar axis (if 
one assumes an isotropic intensity of the produced radiation). 

The apparent luminosity of the accreting star depends on the 
binary inclination i and the orbital phase (f>. We choose <f> = at 
the superior conjunction of the accretor, i.e. when the companion is 
between the observer and the accretor. Then the angular position of 
the observer as viewed from the accretor is given by 9 OOB = i and 

<£obs = — <j>. 

Let CI and i~i bs be unit vectors corresponding to (9, tp) and 
(#obs, fobs), respectively, and fj, = (O • fiobs). The apparent lu- 
minosity is 



(20) 



L obs (i,0)=4 J F(e,<f)fiH(fi)dS 
GMMto 



■kR, 



fiH(fi) dn 



Here H(n) is the Heaviside step function. To compute the integral, 
we substitute n = cos 6 cos i + sin 9 sin i cos(ip + <f>) and take 
9(R t ) and tp(R t ) from equations (12,13). 

Fig. 4 shows L b s found for accretion flows with lo = 0.51, 
and a = 0, 0.5, and 1 (same cases as shown in Fig. 3). One sees 
strong variations of L b s with the orbital phase. The amplitude of 
the variability reaches ~ 3 at i = tt/2 and vanishes at i = 0. 
The maximum at <j> ~ 77r/4 is produced by the bright spot on the 
surface of the star (see Fig. 3). Note that at large i one should also 
take into account the eclipse by the donor. 
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5 CAUSTICS OUTSIDE THE ACCRETOR 
5.1 Formation of caustics 

We now address the case I > U. If the accretion flow is symmetric 
about the equatorial plane then a streamline with / > coming 
from above will collide in this plane with the symmetric streamline 
coming from below. The collision occurs at r — (l/l*) 2 R*. We 
study inflows with A > (see eq. 14), so that there are no inter- 
sections of the ballistic trajectories outside the equatorial plane^. 

The collision in the equatorial plane is associated with a cou- 
ple of shocks that envelope the caustic from above and below. The 



* In general, depending on the distribution £(0oo, foo), such 'early' inter- 
sections are possible. In that case, the increased pressure near A = would 
alter the trajectories. It would cause a relatively modest focusing effect on 
the streamlines, without substantial energy release. By contrast, the even- 
tual collision in the symmetry plane liberates a large fraction of the infall 
kinetic energy. 
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caustic shock is similar to the shock on the surface of a NS (now 
the symmetry plane plays the role of a 'hard surface'). Like the case 
of collision with the star, we assume that the shocks are pinned to 
the caustic. The matter is then in free fall until it reaches the equa- 
torial plane. The flux of mass impinging the caustic determines its 
brightness. 

Note that the occurrence of caustics and the associated energy 
release are not caused by the centrifugal barrier often mentioned in 
the literature. The radial infall is not stopped by rotation when mat- 
ter reaches the caustic. Rather, the two symmetric streams 'miss' 
the center and collide. As a result they cancel 9— components of 
velocity and continue to accrete in the equatorial plane. The subse- 
quent accretion proceeds via a fast disc (see Beloborodov & Illar- 
ionov 2000). Note that the free fall from infinity onto the equa- 
torial plane executes only 1/4 of the full turn around the z-axis 
(ip — ipoo — 7r/2). The centrifugal barrier would stop the infall 
at the periastron radius r p — (l/U) 2 (R*/2) after 1/2 of the full 
turn (see eqs. 11, 13, and Fig. 2). The collision thus happens before 
the centrifugal barrier. 

In the most general case, the streams from above and below 
are not symmetric about the equatorial plane. The caustic may then 
be a time-dependent warped surface. This situation would be the 
subject of a separate study. In this paper, we restrict our consid- 
eration to the inflows which are symmetric about the equatorial 
plane, but not necessary axisymmetric. Then the general shape of 
the caustic is an asymmetric ring in the equatorial plane. 

5.2 The caustic shape 

We use the polar coordinates (r, <p) on the equatorial plane. A 
streamline starting at 9^ , ipoo reaches the plane at 

_ l 2 (0 o 



•p = <Poo + n/2. 



(21) 



GM 

We thus have a mapping (6^ ,fao) — * (r,ip). With a homogeneous 
accretion rate through Soo , we get the flux of matter impinging the 
caustic on one side at given r, ip, 



dM. M tot 



dr 



dcos( 



(22) 



The equatorial streamlines have the highest angular momentum 
^max(Voo) and the outer edge of the caustic is defined by the 
streamlines with #oo — > tt/2, 



ro(<p) = 



/- 



i{<fo 



Poo = ¥ — vr/2. 



(23) 



GM 

If l(0oa,ipoa) is a differentiable function at 6ca — tt/2 then 
dl/dcos 9oo — and dM/dS — > oo at r — > ro(ip), i.e. the flux of 
matter diverges at the outer edge of the caustic. 
As an illustration take the inflow (18). Then 

M tr . / r \ -1/2 



dM 

ds (r '^ 



87rrro 



V r 



R* <r < r (<p) (24) 



where r (<p) = (/g/GM)(l — a sin p) 2 . Note that the r-integrated 
distribution of the accretion rate is dA//d<^> w A/tot /4-7T = 
const at ro 2> R*- At a given ip, the minimum dM/dS = 
(3V3/8)(M to t/27rrg( ¥ j)) is achieved at r min = (2/3)r Q (<p). In 
Fig. 5 we take lo = v3i» and compare the distributions dM/dS 
in the cases of a = and a = 0.5. 

Note two general features of the caustic: 

(i) The asymmetry of 1(8 ^ , fx) in the y-direction results in the 
asymmetry of the caustic in the a>direction. This is a consequence 



(a) 





Figure 5. The distribution of d M/dS over the caustic surface (at r > R* ) 
and the surface of the accretor (at r < R t ). The view is taken from the top, 
i.e., the plane of the figure corresponds to the equatorial plane, 9 = it/2. 
The dashed curve shows the radius i?» . The solid curves display the con- 
tours of constant < : A / /dS with a logarithmic step 0.3. Dark regions corre- 
spond to high dAf /dS. The assumed angular momentum of the accretion 
flow is given by equation (18) with lo = V3U . Plot (a) shows the axisym- 
metric case (a = 0) and plot (b) shows the case a = 0.5. The outer edge of 
the caustic is rn((p) = 3-R*(l — asimp) 2 . 



of the fact that the streamlines execute a 1/4 turn by the moment 
of collision. The caustic thus appears in the front of the accretor 
orbiting the donor star. If the trapped l z had the opposite sign, the 
caustic would appear in the rear of the moving accretor. 

(ii) The outer edge of the caustic is formed by the nearly equa- 
torial streamlines, 9oo — > tt/2. The edge is sharp, with dM/d5 — > 
oo if the /-distribution is smooth (differentiable) at 9^ = tt/2. 
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6 DISCUSSION 

The regime of accretion studied in this paper applies to weakly 
magnetised accretors. In the case of accretion onto a strongly mag- 
netised neutron star, the effective radius of the accretor is the 
Alfvenic radius, Ra and therefore the characteristic h in that prob- 
lem is ~ (GMR A ) 1/2 . 

In reality one expects wind-fed accretion flows to be time- 
dependent. Numerical simulations of the Compton heated subsonic 
region at R <; Rc show that the flow is unsteady (Igumenshchev 
et al. 1993). Fluctuations at R <; Rc imply variable boundary con- 
ditions for the free fall inside Rc- The typical time-scale of the 
variations is of order of the free-fall time at Rc ■ The accretion flow 
is thus expected to fluctuate on time-scales ~ 10 s. 

Throughout the paper we assumed that the flow is symmetric 
with respect to the equatorial (binary) plane. If the symmetry is bro- 
ken, a warped caustic may form near the accretor. The changes in 
the caustic shape are then driven by the ram pressure of the collid- 
ing gas. The warped caustic is likely to be unstable, leading to os- 
cillations/fluctuations on time-scales ~ (7?»/GM) 1/2 ~ 0.1 ms. 

The pattern of accretion studied in this paper assumes that 
the shocks on the star surface/caustics are radiatively efficient and 
pinned to the star surface and/or the equatorial plane. At a low ac- 
cretion rate (which implies low density) the protons heated in the 
shock may find it difficult to pass their energy to the electrons on 
the free-fall time-scale. Then the shocked gas cannot radiate the 
heat and a variable pressure-driven outflow is likely to form. 

On the observational side, the possible modulation of X-ray 
emission with the binary period is especially interesting. Note that 
orbital modulations of soft X-rays are known to occur in wind-fed 
systems as a result of photoelectric absorption in the wind (e.g. Wen 
et al. 1999; Balucihska-Church et al. 2000). The study of orbital 
modulations in the hard X-ray band would be especially helpful 
since they can be caused by intrinsic anisotropy of the source. 
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